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Present photoionization experiments cannot measure molecular frame photoelectron angular dis- 
tributions (MFPAD) from the outermost valence electrons of molecules. We show that details of 
the MFPAD can be retrieved with high-order harmonics generated by infrared lasers from aligned 
molecules. Using accurately calculated photoionization transition dipole moments for fixed-in-space 
molecules, we show that the dependence of the magnitude and phase of the high-order harmonics 
on the alignment angle of the molecules observed in recent experiments can be quantitatively re- 
produced. This result provides the needed theoretical basis for ultrafast dynamic chemical imaging 
using infrared laser pulses. 

PACS numbers: 33.80.Eh, 42.65.Ky 



Photoionization (PI) is the basic process that allows 
the most direct investigation of molecular structure. 
Measurements of total, partial and differential PI cross 
sections have a long history. Until recently, however, 
almost all the experiments are performed from an en- 
semble of randomly oriented molecules. Thus the rich 
dynamical structure of photoelectron angular distribu- 
tion for fixed-in-space molecules predicted in the seminal 
paper by Dill [l| more than 30 years ago still remains 
largely unexplored. In recent years, fixed-in-space PI has 
been investigated with x-ray or VUV photons by using 
the photoelectron-photoion coincidence technique 0, Q • 
Following inner-shell or inner- valence-shell ionization the 
molecular ion dissociates. The molecular axis at the time 
of ionization is inferred from the direction of motion of 
the fragment ion if the dissociation time is short com- 
pared to the rotational period. Clearly this method is 
not applicable to PI from the highest occupied molecular 
orbital (HOMO). 

Consider photoionization of a linear molecule. For a 
polarized light, the differential cross section can be ex- 
pressed in the general form Q (atomic units are used 
throughout, unless otherwise indicated) 
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In this expression, the molecular axis is fixed, the direc- 
tions of light polarization and electron emission are given 
by f2ri and Qk, respectively. This is called the molecular 
frame photoionization angular distribution (MFPAD). 
We calculate the continuum wavefunction ^Jfe/m an( ^ tne 

dipole Ii mil = fc 1/2 (*i| r Ml*/,fcim) ( witn r M = z for lin " 
ear polarization) using the iterative Schwinger variational 
method Q within the complete-active-space configura- 
tion interaction scheme [5|. To compare with measure- 
ments, integration of Eq. (1) over unobserved variables 



has to be performed, thus losing much valuable informa- 
tion on the structure of the molecule. 

In this Letter we show that recent advances in the gen- 
eration of high-order harmonics from molecules by in- 
tense infrared lasers has made it possible to probe fixed- 
in-space molecular PI from HOMO. Gas phase molecules 
can be impulsively aligned by a short sub-picosecond in- 
frared laser [g]. After the pulse is turned off, molecules 
will be partially aligned or anti-aligned at the time in- 
tervals of rotational revival @, Q • During these revivals 
which last for tens or hundreds of femtoseconds, another 
probe pulse can be used to illuminate molecules to ob- 
serve the emission of high-order harmonics. The first 
experiment of this type was reported by Kanai et al. Q 
using N2, O2 and CO2 molecules. Since then many more 

experiments have been reported 00,(3 EI Kill- To 

study the dependence of high-order harmonics generation 
(HHG) on the alignment of molecules, the probe laser can 
be applied at the different time as the molecules evolve, or 
by changing the direction of the probe laser polarization 
at a fixed time delay. More recently, HHG from mixed 
gases 13, II I and interferometry technique [l2j have also 
been used such that the phase of the HHG can also be 
measured. 

HHG is a highly nonlinear process and can be under- 
stood using the three-step model Electrons are 
first released from the molecule by tunnel ionization 
and thrown into the laser field. As the laser's electric 
field changes direction, the electrons may be driven back 
to the ion core after gaining kinetic energy from the 
field. When the electrons recombine with the ion, high- 
order harmonic photons are emitted. Since the photo- 
recombination (PR) in the last step is the time-reversed 
process of PI, it is thus expected that HHG spectra hold 
valuable information on the molecule. Indeed, HHG 
spectra had been used to extract the HOMO in N 2 us- 
ing tomographic method (lEj . However the tomographic 
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method relies exclusively on approximating the contin- 
uum electrons by plane waves which is known to fail in 
general from the PI studies. In we have developed 

a quantitative rescattering theory (QRS) for HHG where 
complex HHG dipole moment D(lo,9) from a fixed-in- 
space molecule can be expressed as a product of the tran- 
sition dipole d{u,6) = E; mM hm^Y* m {n k )Y^{Q, n ) (with 
k || n) and the returning electron wave-packet W(Ek, 9), 



D(u,0) = W{E k ,d) x d(u,9). 



(2) 



This model establishes the relation between the induced 
dipole moment with the transition dipole for PR (or PI) . 
Here 9 is the alignment angle of the molecular axis with 
respect to the polarization direction of the probe laser, 
and u!=I p + Ek, with to the photon energy of the harmon- 
ics, I p the ionization potential, and Ek — fc 2 /2 the "in- 
cident" energy of the returning electron. The validity of 
this equation has been established for atomic targets (no 
9 dependence for atoms) 1(| 17 1 as well as for H J target 
[HI]. In both cases, HHG spectra can be calculated accu- 
rately by solving the time-dependent Schrodinger equa- 
tion (TDSE). In this Letter, we obtained transition dipole 
d(oj, 9) from the state-of-the-art PI calculations. The 
wave-packet W(Ek,9)js obtained from the strong-field 
approximation (SFA)jH, H, [l9[. We then use Eq. (2) 
to obtain D(uj,9) for each fixed-in-space molecule. To 
compare with experiments, proper coherent convolution 
with the partial alignment of molecules has to be carried 
out [see Eq. (2) of Ref. [l(J]. 

In Fig. 1(a), we show the theoretical CO2 PI differen- 
tial cross section for the emission direction along the light 
polarization (fe || n), which are to be compared to HHG 
data below. By changing the orientation of the fixed-in- 
space molecule, complementary information on the an- 
gular distribution of photoelectrons from fixed-in-space 
molecules is obtained. In order to compare with HHG 
data, the photon energies are expressed in units of pho- 
ton energy of 800-nm laser (1.55 eV). First, we note that 
the differential cross section is large when the molecule 
is aligned at a large angle. The cross sections vanish at 
9 = and 7r/2 due to the ir g symmetry of the HOMO and 
the dipole selection rule for the final state. Second, the 
cross section shows near-zero minima for angles around 
35° to 45° for harmonics above the 25th order, or H25. 
The position of the minimum moves to a larger angle, as 
energy increases. These minima have been interpreted as 
evidence of interference of electron waves from the atomic 
centers. In Fig. 1(b), we show the phase of the PI dipole 
vs photon energy, for alignment angles from 20° to 60°. 
Between 30° and 45° the dipole phase undergoes rapid 
change between H25 and H39. This is the region where 
the cross sections are near the minima [see Fig. 1 (a)] . The 
phase change is about 2.0 to 2.5, instead of ir as should be 
expected if the cross section indeed goes to zero. For an- 
gles below 20° and above 60° the phase evolves smoothly 
vs photon energy. The "rich" structures in the theoreti- 
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FIG. 1: (Color online) CO2 differential photoionization cross 
section (a) and phase of the transition dipole (b). The photon 
energies are expressed in units of harmonic orders for 800-nm 
laser, (c) HHG yield as function of angle between pump and 
probe laser polarization directions. The laser intensity and 
duration are of 0.55 x I (I = 10 14 W/cm 2 ), and 120 fs for 
the pump pulse, and 2.5 x Jo and 25 fs for the probe pulse. 
Rotational temperature is taken to be 105 K. (d) Harmonic 
phase (relative to the phase from Kr) for parallel aligned and 
perpendicular aligned ensembles, under two sets of parame- 
ters (solid and dashed lines) that lead to two different degrees 
of alignment distributions. Experimental data from Boutu el 
al. [131 are shown as symbols. 



cal amplitude and phase of d(ui,9) shown in Fig. 1 have 
never been observed directly in PI experiments. Below 
we show that they have been observed in recent HHG 
measurements. 

In Fig. 1(c) we show the simulated typical HHG yields 
from aligned CO2, as a function of the angle between the 
pump and probe lasers polarization directions. These re- 
sults are consistent with recent experiments S 2(J 21 1, 
and with earlier theoretical results 2^, IU 24|. The re- 
sults also resemble the data for the induced dipole re- 
trieved from mixed gases experiments by Wagner et al. 
[HI (see their Fig. 4). In our simulation, the laser param- 
eters are taken from the experiment of fl2| . The align- 
ment distribution is obtained from numerical solution of 
the TDSE within the rotor model 0, 0, [23j| • At the de- 
lay time corresponding to the maximum alignment near 
half-revival (1/2 of the rotational period), a probe beam 
with polarization direction varying from to 90° is used 
to generate high-order harmonics. Comparing Fig. 1(c) 
with Fig. 1(a), we note that the HHG yields follow the 
general angular and photon energy dependence of the 
the differential PI cross sections. Both figures show large 
yields at large angles, and minima near 35° for harmonic 
orders above H31. Due to the averaging over the molec- 
ular alignment distributions, the angular dependence of 
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HHG is smoother. 

We next show that the phase of hxed-in-space PI tran- 
sition dipole can also be probed by comparing to the 
phase of the harmonics. Experimentally the harmonics 
phase can be extracted from measurements of HHG us- 
ing mixed gases 11, 13 1 or interferometry method In 
Fig. 1(d), we show the recent experimental data of Boutu 
et al. [13J where the harmonics phases (relative to that 
from Kr) are obtained for the parallel aligned and per- 
pendicularly aligned ensembles, shown by black and red 
symbols, respectively. For the latter, the phase does not 
change much within H17 to H31. For the parallel aligned 
molecules, the phase jump from H17 to H31 was reported 
to be 2.0 ± 0.6 radians (in the figure it was shown at 2.6 
radians). Our simulation results are shown for two dif- 
ferent ensembles, with alignment distributions confined 
in a cone angle of 25° and 35° at half maximum, re- 
sulting from different pump beams and gas temperatures 
(see Caption). From Fig. 1(b), we note that for a fixed 
harmonic, the phase from the parallel aligned ensemble 
is larger if the angular spread of the molecular distri- 
butions is smaller. Therefore, the position of the phase 
jump for the less aligned ensemble (dashed line) slightly 
shifted toward higher energies. For the perpendicularly 
aligned molecules, the phase is small and does not change 
much with the harmonic order. This is consistent with 
the phase shown in Fig. 1(b). Thus the phase of HHG 
in Fig. 1(d) is seen to mimic the phase of the PI dipole 
shown in Fig. 1(b). 

Next we show in Fig. 2 how the HHG yields change 
with the time delay near the 3/4-revival where the 
molecule can be most strongly aligned. Our simulations 
were carried out with the laser parameters from Zhou 
et al. (l2| . The degree of alignment (top panel), as 
measured by (cos 2 9) is maximum when molecules are 
maximally aligned (vertically in the figure) and minimum 
when molecules are anti-aligned. Our results are shown 



for the same harmonics that had been analyzed in 12 1 
(see their Fig. 2 - according to a new assignment by the 
authors of [12j |. the harmonic orders should be properly 
shifted down by two orders, compared to the ones given 
in their original paper [2l[). For the lower harmonics 
H21 and H25, the yields follow the inverse of the align- 
ment parameter. This is easily understood from Fig. 1(a) 
which shows that the PI cross sections at large angles are 
much larger than at small angles. For the higher orders, 
Fig. 1(a) indicates that the PI cross sections show two 
humps, with the one at smaller angles only a factor of 
about two to three times smaller than the other. Qual- 
itatively, this explains why the HHG yields for H31 and 
up show a pronounced peak for the parallel alignment. 
We note a quantitatively good agreement between our 
calculations and the experimental data. 

Evolution of the HHG yields has also been studied by 
Wagner et al. [10| using the mixed gases technique. The 
yield and phase for H31 vs time delay near half-revival 
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FIG. 2: (Color online) Normalized HHG yield (vs isotropically 
distributed molecules) from CO2 for different harmonics as 
function of pump-probe delay time near 3/4-revival. Laser 
parameters are the same as in Fig. 1(c). The experimental 
data are taken from Zhou et al [l^. The alignment parameter 
(cos 2 6) is plotted for reference (a) . 



are shown in Fig. 3(c) and (b), respectively. Our sim- 
ulations (solid lines) are also in good agreement with 
these measurements. The HHG yield is small when the 
molecules are maximally aligned. This is in agreement 
with H31 seen in Fig. 2(c) for the delay time near 3/4- 
revival. However, we note that the phase is maximum 
when the molecules are maximally aligned. This can be 
understood from the phase in Fig. 1(b) where it shows 
that the phase is large when the alignment angle is small. 
It also shows that the phase is small when the alignment 
angle is large, thus for anti-aligned molecules the phase 
of the HHG should be small, as seen in the experimental 
data and in the simulation. 

We remark that the laser intensity used in the JILA 
experiments is quite high, thus the effect of the depletion 
of the grou nd state by the laser should be considered 
22l 23l 25l ] . For CO2, the alignment dependence of the 
ionization rate is still not fully settled yet. The results 
from the MO-ADK theory the SFA theory and the 
experiment 27J all disagree with each other, with the 
SFA predicting a peak near 40°, compared to 30° from 
MO-ADK theory, and 45° from experiment 27j. Using 
ionization rates from MO-ADK, we have not been able 
to simulate experimental data accurately. We thus used 
ionization rates from SFA, except that we renormalize the 
SFA rate to that of the MO-ADK rate at laser intensity 
of Iq, which gives a factor of 10. Note that the same 
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FIG. 3: (Color online) Alignment parameter (cos 2 9) (a), 
phase (b) and magnitude (c) of the 31th order harmonic from 
CO2, as functions of pump- probe delay time near half-revival. 
The experimental data are taken from Wagner et al [l(J. The 
laser intensity and duration are of 0.38 x Io and 140 fs for 
the pump, and 2.5 x Io and 30 fs for the probe. Rotational 
temperature is taken to be 70 K. 



correction factor has been found for the SFA ionization 
from Kr, which has almost the same ionization potential 
as for C02- With the corrected SFA rate, we found that 
our simulations giv e a better quantitative agreement with 



the JILA data [12 1 



We have also studied N2 molecules. The fixed-in-space 
differential PI cross sections are shown in Fig. 4(a). For 
HI 7 to H21 the cross sections are quite large since this 
is in the tail region of the famous shape resonance in 
N2. For these orders, the cross sections are highly for- 
wardly peaked. For energies above H25, the cross sec- 
tions become much smaller, and they are of the same 
order of magnitude between parallel and perpendicular 
alignments. However, note that the cross sections have 
minima at large angles around 60° - 70°. The (normal- 
ized) HHG spectra of N2 has been measured vs the time 
delay in 0]. Their results are shown in Fig. 4(c), along 
with the results from our simulation. Clearly for H21- 
H25, the yield (normalized to the isotropic distribution 
case) is much larger than that from the higher harmon- 
ics when the molecules are maximally aligned, which is 
consistent with the PI cross sections shown in Fig. 4(a). 
Once again, the experimental HHG spectra can be re- 
produced within the QRS based on the accurate dipole 
transition amplitudes from PI calculation, thus allowing 
us to probe the MFPAD for HOMO. 

In conclusions, we have shown that molecular frame 
photoionization can be probed directly using laser- 
generated high-order harmonics from aligned molecules. 
This method is particularly useful for probing fixed-in- 
space molecular photoionization from the HOMOs. Al- 
ternatively, using accurate PI dipole matrix elements cal- 
culated from the state-of-the-art molecular photoioniza- 
tion codes, we have illustrated that the nonlinear HHG 
spectra can be accurately calculated [see Eq. (2)]. By 



FIG. 4: (Color online) (a) Same as in Fig. 1(a), but for N 2 . (b) 
and (c): Alignment parameter (cos 2 6) and harmonic yields 
H21-H25 and H35-H37 from N2. The experimetal data and 
laser parameters are taken from Itatani et al Q. 



varying laser intensity, HHG spectra from lower-lying or- 
bitals can also be similarly probed [3]. Having estab- 
lished the connection between HHG and the PI dipole 
matrix elements, see Eq. (2), time- resolved PI cross sec- 
tions can be extracted from HHG spectra in a typical 
pump-probe experimental setup, thus opening up the op- 
portunity of using infrared laser pulses for ultrafast dy- 
namic chemical imaging. 
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